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I. INTRODUCTION 

Atomic helium and light helium-like ions have long 
been attractive subjects of theoretical and experimen- 
tal investigations. From the theoretical point of view, 
helium-like atoms are the simplest few-body systems. As 
such, they are traditionally used as a testing ground for 
different methods of the description of atomic structure. 
On the experimental side, small natural linewidths of 
transitions between the metastable 3 P and 3 S states of 
helium-like ions permit spectroscopic measurements of 
high precision. For atomic helium, experimental investi- 
gations are nowadays carried out with the relative accu- 
racy up to 7 x 1CP 12 [1]. An advantage of the helium- 
like ions as compared to, e.g., hydrogen-like ones, is that 
the transition frequency increases slowly with the nuclear 
charge number Z (~ Z). This feature ensures that wave- 
lenghts of a significant part of the helium isoelectronic 
sequence fall in the region suitable for accurate experi- 
mental determination. 

There are presently two main theoretical approaches 
that allow one to systematically account for the electron- 
correlation, relativistic, and quantum electrodynamical 
(QED) effects in few-electron systems. The first one, 
traditionally used for light systems, relies on an expan- 
sion of the relativistic and QED effects in terms of a 
and Za (a is the fine-structure constant) and treats the 
nonrelativistic electron-electron interaction nonperturba- 
tively. This approach started with the pioneering works 
of Araki [2] and Sucher [3] , who derived the expression 
for the Lamb shift in many-electron systems complete 
through the order ma 5 . The other approach aims pri- 
marily at high-Z ions. It docs not use any expansion 
in the binding-strength parameter Za (and thus is of- 
ten referred to as the all-order approach) but treats the 
electron-electron interaction within the perturbative ex- 
pansion with the parameter X/Z. A systematic formula- 
tion of this method is presented in Ref. [4] . 

These two approaches can be considered as comple- 
mentary, the first being clearly preferable for light atoms 
and the second, for heavy ions. The intermediate region 



of nuclear charges around Z = 12 is the most difficult 
one for theory, as contributions not (yet) accounted by 
either of these methods have their maximal value there. 
In order to provide accurate predictions for the whole 
isoelectronic sequence, it is necessary to combine these 
two approaches. 

For the first time a combination of the complementary 
approaches was made by Drake [5]. His results for en- 
ergies of helium-like ions comprise all effects up to order 
ma 5 in the \ow-Z region, whereas in the high-Z region, 
they are complete up to the next-to-the-leading order 
in 1/Z for nonradiative effects and to the leading or- 
der, for radiative effects. Since then, significant progress 
was achieved in theoretical understanding of energy lev- 
els of atomic helium, whose description is now complete 
through order ma 6 [6, 7]. Also in the high-Z region, 
theoretical energies have recently been significantly im- 
proved by a rigorous treatment of the two-electron QED 
corrections [8], which completed the 0(1/Z) part of the 
radiative effects. 

In the present investigation we aim to improve theo- 
retical predictions of the n — \ and n — 2 energy levels of 
light helium-like ions. To this end, we perform a calcula- 
tion that includes all QED and recoil effects up to orders 
ma 6 and m 2 /M a 5 (M is the nuclear mass). In order to 
establish a basis for merging the current approach with 
the all-order calculations, we perform an extensive anal- 
ysis of the 1/Z expansion of individual corrections. This 
analysis also provides an effective test of consistency of 
our calculational results and of the 1/Z-expansion data 
available in the literature. 



II. THEORY OF THE ENERGY LEVELS 

In this section, we present a summary of contributions 
to the energy levels of two-electron atoms complete up to 
orders ma 6 and m 2 /M a 5 . 

According to QED theory, energy levels of atoms are 
represented by an expansion in powers of a of the form 



E(a) = £?« + E^ + E^ + + E W 



(1) 
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where = ma n £^ is a contribution of order a n and 
may include powers of In a. Each of £^ is in turn ex- 
panded in powers of the electron-to-nucleus mass ratio 
m/M 
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where £ M denotes the correction of first order in m/M 

and £^1 is the second-order correction. To note, for the 
nonrelativistic energy, it is more natural to expand in 
m r /M (where m r is the reduced mass) rather than in 
m/M, since such expansion has smaller coefficients. For 
the relativistic corrections, however, the natural recoil 
expansion parameter is m/M, so for consistency we use 
it for the nonrelativistic energy as well. 

The terms of the double perturbation expansion (1) 
and (2) are expressed as expectation values of some effec- 
tive Hamiltonians (in some cases, of nonlocal operators) 
and as second- and higher-order perturbation corrections 
induced by these Hamiltonians (operators). It is note- 
worthy that the expansion (1) is employed also for the 
states that are mixed by the relativistic effects, namely 
2 1 Pi and 2 3 Pi. The mixing effects are treated pertur- 
batively. (So, the leading effect due to the 2 1 P\ — 2 3 Pi 
mixing appears naturally as the second-order ma 6 cor- 
rection, together with contributions from other interme- 
diate states.) This differs from the approach used, e.g., 
in Ref. [5], where a two-by-two matrix was constructed 
for this pair of states and the energies were obtained by 
a diagonalization. 

(2) 

The leading contribution to the energy £^o' = £q is the 
eigenvalue of the nonrelativistic Hamiltonian, 



rj + ^r ab 
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The first- and second-order recoil corrections to the non- 
relativistic energy are given by 
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where 



(5) 
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is the mass polarization operator. 



The leading relativistic correction £^ is given by the 
expectation value of the Breit-Pauli Hamiltonian 



[9], 
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The finite nuclear mass correction to the Breit contri- 
bution is conveniently separated into the mass scal- 
ing, the mass polarization, and the operator parts. The 
mass scaling prcfactor is (m r /m) for the first term in 
Eq. (7) and (m r /m) 3 , for all the others. The mass po- 
larization part represents the first-order perturbation of 

£00 by the mass-polarization operator (6). The opera- 
tor part is given by the expectation value of the recoil 
addition to the Breit-Pauli Hamiltonian, 
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£co is the leading QED correction [2, 3]. We divide it 

(5) 

into the logarithmic and the nonlogarithmic parts, £^, = 
£c?(log) + £^(nlog), which are given by 

£<o 5) (log) =^ln(Za)^> 3 (r o6 )) 
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where 



Qab = ^^r+S 3 (r ab ) \nZ 
and the singular operator r -3 is defined by 
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where 7 is the Euler constant. The Bethe logarithm is 
defined as 



ln(fco) 
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The operator Hf s is the anomalous magnetic moment 
correction to the spin-dependent part of the Brcit-Pauli 

(5) 

Hamiltonian. Hf s does not contribute to the energies 
of the singlet states and to the spin-orbit averaged levels 
but it yields the ma 5 contribution to the fine structure 
splitting. It is given by 
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We note that despite the presence of terms with In Z in 
Eq. (10), the correction £<^(nlog) does not have loga- 
rithmic terms in its 1/Z expansion. 

The recoil correction consists of four parts [10], 



4»^(A+^^)l(<J, (15) 



where 
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with Q a defined analogously to Eq. (11), and (m/M)£ 3 

(5) 

is the first-order perturbation of £oo due to the mass- 
polarization operator (6). The operator Hfp ICC yields a 



nonvanishing contribution to the fine-structure splitting 
only. It is given by 
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We note that the last term in Eq. (16) was omitted in 
the original derivation of Ref. [10]. 

The complete result for the m a 6 correction £ ^ to the 
energy levels was derived by one of the authors (K.P.) in 
a series of papers [6, 7, 11, 12] 
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The first term in the above expression contains the com- 
plete logarithmic dependence of the ma 6 correction. 
The part of it proportional to In a was first obtained in 
Ref. [13]. The remaining logarithmic part proportional 
to InZ was implicitly present in formulas reported in 
Ref. [6, 7] (it originates from the expectation value of the 
operator l/rjjj). In Eq. (19), we group all logarithmic 
terms together so that the remaining part does not have 
any logarithms in its 1/Z expansion. 

The term E scc in Eq. (19) is the second-order perturba- 
tion correction induced by the Breit-Pauli Hamiltonian. 
(More specifically, it is the finite residual after separat- 
ing divergent contributions that cancel out in the sum 
with the expectation value of the effective ma 6 Hamilto- 
nian.) The first part of the effective Hamiltonian, H^) ad , 
originates from the non-radiative part of the electron- 
nucleus and the electron-electron interaction. The next 
two terms, and H^, are due to the one-loop and 

two-loop radiative effects, respectively. The last two 
parts and H^ mm are the spin-dependent opera- 

tors first derived by Douglas and Kroll [14]. They do not 
contribute to the energies of the singlet states and to the 
spin-orbit averaged levels. Expressions for these opera- 
tors are well known and are given, e.g., by Eqs. (3) and 
(7) of Ref. [15]. The non-radiative part of the ma 6 effec- 
tive Hamiltonian is rather complicated. For simplicity, 
we present it specifically for a two-electron atom. The 
corresponding expression reads [6, 7] 
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where P = pi + p 2 , p = (pi - p 2 )/2, f = f 1 - f 2 . We 

note that the operator P^ad i s defined in such a way that 
its expectation values does not contain any logarithmic 
terms in the 1/Z expansion, as the last term of Eq. (20) 
is compensated by the corresponding contribution from 
the 1/r 3 operator. 

The effective Hamiltonians induced by the radiative 
effects are [6, 16, 17] 
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The second-order correction can be represented as 
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where P^ and are the spin-independent and spin- 
dependent parts of the Breit-Pauli Hamiltonian (7), re- 
spectively. The operator is obtained from by a 
transformation that eliminates divergences in the second- 
order matrix elements [6]. The transformed operator is 
given by 



5^ + 
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where V 2 V 2 . is understood as a differentiation of the 
wave function on the right hand side as a function (omit- 
ting 8 3 (r)) and V = -Z/n - Z/r 2 + 1/r. 
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An intriguing feature of the formulas presented in this 
section is that the logarithmic dependence of all of them 
enters only in the form of ln(Za). This is not at all obvi- 
ous a priori since ln(Za) appears naturally only in contri- 
butions induced by the electron-nucleus interaction. The 
effects of the electron-electron interaction usually yield 
logarithms of a, whereas logarithms of Z are implicitly 
present in matrix elements of singular operators. The 
fact that logarithms of a and logarithms of Z have the 
coefficients that match each other comes "accidentally" 
from the derivation. 

The complete result for the corrections of order m a 7 
for the helium Lamb shift is not presently available (it 
is known for the fine-structure splitting only [15, 18]). 
One can, however, easily generalize some of the hydro- 
genic results, namely those that are proportional to the 
electron density at the nucleus. These are (i) the one- 
loop radiative correction of order ma {Zaf \n 2 (Za)- 2 , 
(ii) the two-loop radiative correction of order ma 2 (Za) 5 , 
and (iii) the nonrelativistic correction due to the finite 
nuclear size. The first two effects yield the main contri- 
bution to the higher-order remainder function of S states 
in light hydrogen-like ions. We expect that they domi- 
nate for light helium-like ions as well. 

Following Ref . [5] , we approximate the higher-order ra- 
diative ("rad") and the finite-nuclear-size ("fs") correc- 
tion to the energies of helium-like ions by 



c(7+) = -(7+) (Ej S 3 (n)) 



(25) 
(26) 



where the subscript H corresponds to the "hydrogenic" 
limit, i.e., the limit of the non- interacting electrons, and 



<E«vo>* = v( 1+ 



(27) 



The approximation of Eqs. (25) and (26) is exact for the 
above-mentioned corrections proportional to the electron 
density at the nucleus. It is expected also to provide 
a meaningful estimate for contributions that weakly de- 
pend on n (such as the nonlogarithmic radiative correc- 
tion of order ma(Za) 6 ). Moreover, this approximation 
is exact to the leading order in the 1/Z expansion, thus 
providing a meaningful estimate for high-Z helium-like 
ions as well. 

For all the states under consideration except 2 1 P 1 and 
2 3 Pi, the "hydrogenic" remainder function is just the 
sum of the corresponding remainders for the two elec- 
trons in the configuration, 



c(7+) 



4 7 d +) (i*) 



■4 7 d + Vi). 



(28) 



For the 2 1 P\ and 2 3 Pi states, the Dirac levels need to 
be first transformed from the jj to the LS coupling and 



thus [5] 



Ct^Pi) = ^(1*) + 3 ^(2p 3/2 ) + §^ ; (2pi/ 2 ) • 

(29) 

In our calculation, the one-electron remainder func- 
tion £ rad (nlj) includes all known contributions of order 
ma 7 and higher coming from (i) the one-loop radiative 
correction, (ii) the two-loop radiative correction, (iii) the 
three-loop radiative correction, see the review in Ref. [19] 
and Ref. [20] for an update on the two-loop remainder 
function. 

Besides the finite-nuclear-size and radiative correc- 
tions, there are also non-radiative effects, denoted as 
4ad an d estimated within the 1/Z expansion. More 
specifically, we include the higher-order remainder due 
to the one-electron Dirac energy and due to the one- 
photon exchange correction. They enter at the order 
ma 8 only but are enhanced by factors of Z 8 and Z 7 , 
respectively. Despite this enhancement, numerical con- 
tributions of these effects are rather small for the ions 
considered in the present work. 



'(7+) 



(7+), 



(7+) 
rad 
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III. RESULTS AND DISCUSSION 

A. Numerical results 

The nonrelativistic energies and wave functions are ob- 
tained by minimizing the energy functional with the ba- 
sis set constructed with the fully correlated exponential 
functions. The choice of the basis set and the general 
strategy of optimization of the nonlinear parameters fol- 
low the main lines of the approach developed by Ko- 
robov [21, 22]. The calculational scheme is described in 
previous publications [6, 7, 15] and will not be repeated 
here. Numerical values of the nonrelativistic energies of 
helium-like ions with the nuclear charge Z = 3 ... 12 are 
presented in Table I. The results were obtained with 
N = 2000 basis functions and are accurate to about 18 
decimals (more than shown in the table). The energy 
levels of the helium atom traditionally attract special at- 
tention, so we present the corresponding results in full 
length below. Our numerical values of the upper varia- 
tional limit of the nonrelativistic energies of helium are 



^>(2 3 S) = 

egpp) = 

£^(2 3 P) = 



-2.903 724377 034 119 598310 ±£, (30) 

-2.145 974 046 054 417 415 799 ±%, (31) 

i 

(32) 



-2.175 229 378 236 791 305 738 977 ±2 , 



-2.123 843 086 498 101 359 246 ^ , (33) 
-2.133 164 190 779 283 205 146 96±? . 



(34) 



6 



The value for the ground state is given only for com- 
pleteness, since much more accurate numerical results are 
available in the literature [22, 23]. The numerical results 
for the leading relativistic correction £^ are summarized 
in Table II. 

The leading QED correction is given by Eqs. (9), 
(10), and (15). Computationally the most problematic 
part of it is represented by the Bcthe logarithm ln(fc ) 
and its mass-polarization correction ln(fco)M- Accurate 
calculations of ln(fc ) were performed by Drake and Gold- 
man [24] for helium-like like atoms with Z < 6 and by 
Korobov [25] for Z = 2. Calculations of the recoil correc- 
tion to the Bethe logarithm were reported by Pachucki 
and Sapirstein [10] for Z = 2 and by Drake and Gold- 
man [24] for Z < 6. In the present investigation, we per- 
form accurate evaluations of the Bethe logarithm ln(fco) 
and its recoil correction ln(/co)M for helium-like ions with 
Z < 12. The calculational approach is described in Ap- 
pendix A. 

Table III summarizes the numerical results obtained 
and gives a comparison with the previous calculations. 
Numerical values for the Bethe logarithm are presented 
for the difference ln(fco)— ln(Z 2 ) since this difference has a 
weak Z-dependence and does not contain any logarithms 
in its 1 jZ expansion. The table also lists the coefficients 
of the l/Z expansion of ln(fco/Z 2 ). The leading coef- 
ficient c is known from the hydrogen theory; accurate 
numerical values can be found in Ref. [26]. The higher- 
order coefficients were obtained by fitting our numerical 
data. It is interesting to compare them with the analo- 
gous results reported previously by Drake and Goldman 
[24]. For the next-to-the- leading coefficient ci, the re- 
sults agree up to about 4-5 digits for £ states and about 
3-4 digits for P states. For the higher-order coefficients, 
the agreement gradually deteriorates. However, the re- 
sults for the sum of the two expansions agree very well 
with each other. More specifically, the maximal absolute 
deviation between the values of the Bethe logarithms for 
Z > 12 obtained with our 1/Z-expansion coefficients and 
with those by Drake and Goldman is 1 x 10~ 8 for the l 1 ^ 
state, 3 x 10~ 8 for the 2 X S state, 6 x 10~ 9 for the 2 3 S' 
state, 1 x 10~ 7 for the 2 X P state, and 6 x 10~ 8 for the 2 3 P 
state. So, the accuracy of these expansions is sufficient 
for most practical purposes. 

Another part of the calculation of £^ that needs a 
separate discussion is the evaluation of the expectation 
value of the singular operator 1/r 3 , which is defined by 
Eq. (12). The calculational approach is described in Ap- 
pendix B. Total results for the logarithmic and the non- 
logarithmic part of the leading QED correction are sum- 
marized in Tables IV and V, respectively. The results are 
in good agreement with the previous calculations [5]. 

Table VI presents the numerical values of the ma 6 cor- 
rection, the main result of this investigation. The corre- 
sponding calculations for atomic helium were reported 
in Refs. [6, 7]; our present numerical values agree with 
the ones obtained previously. Calculations performed in 
this work for helium-like ions were accomplished along 



the lines described in Refs. [6, 7]. Here we only note that 
calculations for higher values of Z often exhibit a slower 
numerical convergence (and numerical stability) than for 
helium, especially so for the second-order corrections in- 
volving singular operators. The variational optimization 
of nonlinear parameters for the symmetric second-order 
corrections was performed in several steps with increas- 
ing the number of basis functions on each step up to 
N = 1000 or 1200. The final values were obtained with 
merging several basis sets and enlarging the number of 
functions up to N = 5000 — 7000. The nonsymmetric 
second-order corrections were evaluated as described in 
Ref. [15]. The calculations were performed in the quadru- 
ple, sixtuple, and octuple arithmetics implemented in 
Fortran 95 by V. Korobov [27]. 

Table VII presents the results for the finite nuclear size 
correction and approximate values of the higher-order 
(ma 7 and higher) correction to the ionization energy. 
The uncertainty of the total theoretical prediction origi- 
nates from the higher-order radiative effects; it was esti- 
mated by dividing the absolute value of this correction by 
Z. The values for the root-mean-square radius of nuclei 
were taken from Ref. [28] . 

B. Comparison with the all-order approach 

In this subsection we discuss the calculational results 
obtained for the m a e correction in more detail and 
make a comparison with the results obtained previously 
within the all-order, 1/Z-expansion approach. The log- 
arithmic part of the correction, £( 6 )(log), behaves as 
ma 3 (Za) 3 for large Z and thus corresponds to diagrams 
with three photon exchanges that have not yet been ad- 
dressed within the all-order approach. The nonlogarith- 
mic part £( 6 '(nlog), however, contains some pieces that 
are known and identified below. 

For all states except 2 1 Pi and 2 3 Pi, the leading term of 
the l/Z expansion of £ ^ 6 ^(nlog) is of order m (Za) 6 and 
comes from the Za expansion of the one-electron Dirac 
energy. For the 2}P\ and 2 3 P\ states, the leading term is 
of the previous order in l/Z , m (Za) e Z, and is due to the 
mixing of these levels. More specifically, the mixing con- 
tribution is 5E mix = |(2 1 Pi| J ff( 4 )|2 3 Pi)| 2 /[^o(2 1 Pi) - 
E (2 3 P{)} for the 2 1 P\ state and that with the opposite 
sign, for 2 3 Pi. The contribution of order m(Za) 6 for 
the mixing states comes from the expansion of the one- 
electron Dirac energy and from the expansion of <5P m i x . 

The next term of the l/Z expansion is of order 
ma(Za) 5 and comes from the one-electron one-loop ra- 
diative correction and from the one-photon exchange cor- 
rection. The radiative part is well known [19]. The part 
due to the one-photon exchange was obtained for the l 1 ^, 
2 3 S", 2 3 P , and 2 3 P 2 states analytically in Ref. [29] and 
for the other states numerically in this work. For the 
2 1 Pi and 2 3 Pi states, there is a small additional mixing 
contribution, which we were unable to determine unam- 
biguously. 
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The first two coefficients of the 1/Z expansion of 
£( 6 )(nlog) are listed in Table VI. A fit of our numeri- 
cal data agrees well with these coefficients. The agree- 
ment observed shows consistency of our numerical results 
with independent calculations at the level of the one- 
photon effects. We now turn to the contribution of order 
ma 2 (Za) 4 . This contribution is induced by nontrivial 
two-photon effects, so that a comparison drawn for this 
part will yield a much more stringent test of consistency 
of different approaches. 

The part of £( 6 )(nlog) that is of order ma 2 (Za) 4 is 
implicitly present in the two-electron QED contribution 
calculated numerically in Rcf. [8] to all orders in Za. 
This contribution can be represented as (see Eq. (72) of 
Rcf. [8]) 

A£^f D = ma 2 (Za) 3 \a 31 ln(Za)- 2 + a 30 + (Za) G(Z) 

(35) 

where the remainder function G(Z) incorporates all 
higher orders in Za. The two-electron QED correction 
comprises the so-called screened self-energy and vacuum- 
polarization contributions and the part of the two-photon 
exchange correction that is beyond the Breit approxima- 
tion. 

The coefficients 031 and 030 in Eq. (35) correspond to 
the second term of the 1/Z expansion of the leading QED 

correction . More specifically, 031 corresponds to the 
coefficient c\ from Table IV and 030, to that from Ta- 
ble V. The Z = limit of the higher-order remainder 
function G(Z) corresponds to the third coefficient of the 
1/Z expansion of the correction £( 6 )(nlog), G(0) = C2, 
for all states except 2 1 P i and 2 3 Pi. The values of c 2 
obtained by fitting our numerical data are listed in Ta- 
ble VI. For the 2 1 Pi and 2 3 Pi states, the coefficient C2 is 
not directly comparable with the all-order results because 
of the mixing effects. 

The higher-order remainder function G(Z) inferred 
from the numerical results of Ref. [8] is plotted in Fig. I, 
together with its limiting value at Z = obtained by 
a fit of our numerical data. It should be stressed that 
the identification of the remainder implies a great deal 
of numerical cancellations, especially for the all-order re- 
sults. The comparison drawn in Fig. 1 provides a strin- 
gent cross-check of the two complementary approaches. 
The visual agreement between the results is very good 
for the S states, whereas for the P states a slight dis- 
agreement seems to be present. 

It is tempting to merge the all-order and the Za- 
expansion results by fitting the all-order data for G(Z) 
towards lower values of Z. However, we do not attempt 
to do this in the present work. The reasons are, first, 
that the numerical accuracy of the all-order results is 
apparently not high enough and, second, that the ex- 
pansion of the remainder function G(Z) contains terms 
(Za) ln 2 (Za) and (Za) \n(Za), which cannot be reason- 
ably fitted with numerical data available in the high-Z 
region only. 



C. Total energies 

Our total results for the ionization energy of the n = 1 
and n — 2 states of helium-like atoms with the nuclear 
charge Z = 2 ... 12 are listed in Table VIII. The fol- 
lowing values of fundamental constants were employed 
[19], Poo = 10 973 731.568 527(73) m" 1 and a" 1 = 
137.0359 999 679(94). The atomic masses were taken 
from Ref. [30]. 

The results for atomic helium presented in Table VIII 
differ from those reported previously only because of 
the different approximate treatment of the higher-order 
(ma 7 and higher) contribution employed in this work. 
For the S states of helium, the present values are prac- 
tically equivalent to those of Refs. [6, 7]. (The difference 
is that now we include some contributions of order m a 8 
and higher, which are negligible for helium but become 
noticeable for higher- Z ions.) However, for the helium 
P states, our present estimate of the higher-order contri- 
bution is by about I MHz higher than that of Ref. [7]. 
The reason is that the one-electron radiative correction 
of the p electron state was previously not included into 
the approximation (25). It is included now [see Eq. (29)] 
in order to recover the correct asymptotic behaviour of 
the radiative correction in the high-Z limit. 

A selection of our theoretical results for transition en- 
ergies is compared with the theory by Drake [5, 31] and 
with experimental data in Table IX. Agreement between 
theory and experiment is excellent in all cases studied. 
We observe a distinct improvement in theoretical accu- 
racy as compared to the previous results by Drake. This 
improvement is due to the complete treatment of the cor- 
rections of order ma 6 and m 2 /Ma 5 accomplished in this 
work. 

Theoretical results for the fine structure splitting in- 
tervals 2 3 P — 2 3 Pi and 2 3 Pi — 2 3 P 2 are not analyzed 
in the present work. This is because these intervals can 
nowdays be calculated more accurately (complete up to 
order ma 7 ), like it was recently done for helium [15]. We 
intend to perform such a calculation in a subsequent in- 
vestigation. 

Among the results listed in Table VIII for helium-like 
ions, the ground-state energy of the carbon ion is of par- 
ticular importance, because it is used in the procedure 
of the adjustment of fundamental constants [32] for the 
determination of the mass of 12 C 4+ and, consequently, 
of the proton mass from the Penning trap measurement 
by Van Dyck et al. [33] . Our result for the ground-state 
ionization energy of helium-like carbon is 



P( 12 C 4+ ) = -3162 423.60(32) cm" 



(36) 



which is in agreement with the previous result by Drake 
[5] of -3 162 423.34(15) cm" 1 . Wc note that despite the 
fact that our calculation is by an additional order of a 
more complete than that by Drake, our estimate of un- 
certainty is more conservative. 

In summary, significant progress has been achieved 
during last decades both in experimental technique and 
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theoretical calculations of helium-like atoms. In the 
present investigation, we performed a calculation of the 
energy levels of the n = 1 and n = 2 states of light 
helium-like atoms with the nuclear charge Z = 2 ... 12, 
within the approach complete up to orders ma 6 and 
m 2 /Ma 5 . An extensive analysis of the 1/Z expansion 
of individual corrections was carried out and comparison 
with results of the complementary approach was made 
whenever possible. Our general conclusion is that the 
results obtained within the approaches based on the Za 
and the 1/Z expansion are consistent with each other 
up to a high level of precision. However, further im- 
provement of numerical accuracy of the all-order, 1/Z- 
expansion results and their extension into the lower- Z 
region is needed in order to safely merge the two comple- 
mentary approaches. 
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Appendix A: Bethe logarithm 

Following the approach of Refs. [34, 35], the Bethe log- 
arithm (13) is expressed in terms of an integral over the 
momentum of the virtual photon, 



1 



ln(fc ) = — lim 



r K 

(V 2 ) K + D ln(2AT) + / dkkJ{k) 
Jo 

(Al) 

where D = 2TrZ(5 3 (r 1 ) + 5 3 (r 2 )), V = Vi + V 2 , and 

1 



E -H -k 



(A2) 



For the purpose of numerical evaluation, the integration 
over the photon momentum k is divided into two regions 
by introducing the auxiliar parameter k, 



ln(fc ) = R(k) + -!- / dkk J{k) + [ 
D Jo J K 



dk^f, (A3) 



where the function w(k) represents the residual obtained 
from J(k) by removing all known terms of the large-Zc 
asymptotics, 

1,3 h.2 

w{k) = — J(k) + — (W 2 ) + k -2V2 Zk 1 / 2 + 2Z 2 Ink, 



D 



(A4) 



and R(k) is a simple function obtained by integrating out 
the separated asymptotic expansion terms, 



(V 2 ) , ,0 \ 4y/2Z 2Z 2 (\uk + 1) 
#(«) = K ~FT~ + ln ( 2 «) + 



D 



,1/2 



(A5) 



The calculational scheme employed for the evaluation 
of Eq. (A3) is similar to that previously used [15] for the 
relativistic corrections to the Bethe logarithm. At the 
first step, a careful optimization of nonlinear basis-set 
parameters was carried out for a sequence of scales of the 
photon momentum: ki — 10* and i — l,...,i max , with 
j m = 5 for the S states and i max = 4 for the P states. 
The optimization was performed with incrementing the 
size of the basis until the prescribed level of convergence 
is achieved for the function w(k). At the second step, the 
integrations of the photon momentum k were performed. 
For a given value of k, the function J{k) was calculated 
with a basis obtained by merging together the optimized 
bases for the two closest ki points, thus essentially dou- 
bling the number of the basis functions. The function 
w(k) was obtained from J(k) according to Eq. (A4). 

The integral over fee [0, n] was calculated analytically, 
by diagonalizing the Hamiltonian matrix and using the 
spectral representation of the propagator. The value of 
the auxiliar parameter k was set to k = 100. The integral 
over k <G [k, 00) was separated into two parts, k < 10 lmax 
and k > 10* max . The first part was evaluated numeri- 
cally by using the Gauss-Legendre quadratures, after the 
change of variables t = 1/k 2 . The second part was ob- 
tained by integrating the asymptotic expansion of the 
function w(k). The coefficients of this expansion were 
obtained by fitting the numerical data for w(k) to the 
form 



(A6) 



where pol(a;) denotes a polynomial of x. The total num- 
ber of fitting parameters was about 9—11. The range of 
k to be fitted and the exact form of the fitting function 
were optimized so as to yield the best possible results for 
the known asymptotic expansion terms of J(k). 

The first-order perturbation of the Bethe logarithm by 
the mass-polarization operator can be represented as [10] 



ln(fco) 



M 



Rm(k) + 



dk k Ju{k) 



/ 



k 2 



(A7) 



where 
J M {k) = 2 

V 



E -H -k 



1 



E Q -H Q -k 



[SE-pt -p 2 ] 



1 



E -H -k 



(A8) 
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and 5E = (p\ ■ p 2 )- The perturbed wave function 5(f) is 
defined by 



\S4>) = \8m4>) - |0) 



S M D 
D 



(A9) 



where 5m stands for the first-order perturbation induced 
by the operator p\ ■ p 2 . The asymptotic expansion of 
Jm(^) is much simpler than that of J(k) and is 
just 



Correspondingly, the function Rm{k) is 

2k 



Rm{k) 



I) 



|V 2 |^). 



(A10) 



(AH) 



The numerical evaluation of Eq. (A7) was performed in 
a way similar to that for the plain Bethe logarithm. In 
particular, the same sets of optimized nonlinear param- 
eters were used. Since a high accuracy is not needed 
for this correction, a somewhat simplified calculational 
scheme was used in this case. The high-energy part of 
the photon-momentum integral, k € [100, oo), was eval- 
uated by integrating the fitted asymptotic expansion for 
WM{k), which was taken to be of the form (A6) with 6 — 9 
fitting parameters. 



Appendix B: Expectation value of 1/r 3 

The definition of the expectation value of the regular- 
ized operator 1/r 3 is given by Eq. (12). With the basis- 
set representation of the wave function employed in this 
work, a typical singular integral to be calculated is 



L = 



16tt 2 



cxpC-arx- fo-Tr) 



0(r-e). 
(Bl) 



The straightforward way is to evaluate this integral an- 
alytically for a finite value of the regulator e and then 



expand it in small e. This way is possible, but we prefer 
to use a simpler procedure, which is also the closest to 
the method of evaluation of the regular integrals. 

We recall that all regular integrals are immediately ob- 
tained from the master integral 



1 J ^ eM-an - (3r 2 jr) 



16tt 2 



r\r 2 r 



(a + /3)(/3 + 7 )(a + 7 ) 



(B2) 



by formal differentiation or integration with respect to 
the corresponding parameters. The differentiation over 
a and f3 and an integration over 7 delivers a result for 
the integral of the type 1/r 2 . This integral is conver- 
gent, so the result is exact. The second integration over 
7 (which would yield an integral of the type 1/r 3 ) is di- 
vergent. The simplest way to proceed is as follows. We 
introduce a cutoff parameter for large values of 7, eval- 
uate the integral over 7, and drop all cutoff-dependent 
terms. The expression obtained in this way differs from 
the correct one by a 7 -indcpcndent constant only, which 
can be proved by differentiating with respect to 7. 

The missing constant is most easily recovered by con- 
sidering the behavior of the integral I when 7 — > 00. For 
very large 7, only the region of very small r contributes 
and we have 



00 c 00 



I £ = 2 / drr dr\ r\ / dr 2 r 2 



ri+r 



-ari—pr2—~fr 



|n-r| 



dr- 



(a + /3) 3 j, 
Therefore, 



(B3) 



7 rcg = lim [l e + 7 + In e] 7 =°° - - fl^ In 7 . (B4) 



{a + Pf 



The above equation yields the necessary condition for 
determining the missing constant term in the general ex- 
pression for the regularized integral / rcg . 
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(21 (21 (21 (2) 

TABLE I: Nonrelativistic energies of helium- like ions £ oo , £ M , and £ M i- For the helium atom, the nonrelativistic energy S^,' 

/'2\ (2) 

is given in the text, see Eqs. (30)-(34). For £^ and S M , the results of fitting the numerical data to the form £ = ^ i a/Z 1 are 
also presented. The 1/Z expansion of the second-order recoil correction £^{ was not studied since this correction is relevant 
for the light atoms only. Atomic units are used. 

Z ' l x 5 " 2 1 S " 2 3 S 2 X P 2 3 P 



c(2) 


+ Z 2 (l + l/n 2 )/2 










3 


1.720 086 587 330 694 


0.584 123 254 404 560 


0.514 272 627429 260 


0.631648 922 219 983 


0.597 284 318 602 632 


4 


2.344 433 761576 413 


0.815 126 104 651679 


0.702 833 410 222 385 


0.889 228 377083 556 


0.825 026 856 929 027 


5 


2.969 028 419 757 218 


1.046 471967118 562 


0.891 102 651 185 767 


1.147 716 734 692 201 


1.051862 307 786 520 


6 


3.593 753 398 101 470 


1.277 982 298 711029 


1.079 244 097692 043 


1.406 667686611591 


1.278 289 303 511949 


7 


4.218 554 851227 295 


1.509 584 284 500 004 


1.267 318 262 508 963 


1.665 883 599 383 315 


1.504 498 255 070109 


8 


4.843 404 877 242 074 


1.741242 692 371423 


1.455 352 679 914 990 


1.925 264 764124 369 


1.730 577 283 616 668 


9 


5.468 287636 040 509 


1.972 938 360 878 370 


1.643 361670481692 


2.184 755 723 292 205 


1.956 572 706 521290 


10 


6.093193 484 962 451 


2.204 659 953 768 328 


1.831353 415 926 627 


2.444 323 260 421919 


2.182 511 179 215 097 


11 


6.718116 223 927 278 


2.436 400 325 538 931 


2.019 332 929131276 


2.703 946 297095 311 


2.408 409121058 957 


12 


7.343 051687 353 070 


2.668 154 743 908 730 


2.207 303 452 397112 


2.963 610 824 902 483 


2.634 277196 002 167 


1/Z 


expansion coefficients 










C-l 


5/8 


169/729 


137/729 


1705/6561 


1481/6561 


C() 


-0.157 666 43 


-0.114510 14 


-0.047 409 30 


-0.157 028 66 


-0.072 998 98 


Cl 


0.008 699 03 


0.009 327 61 


-0.004 872 28 


0.026 106 26 


-0.016 585 30 


C2 


-0.000 888 69 


-0.001284 99 


-0.003 457 75 


0.005 782 46 


-0.010 353 67 


C3 


-0.001 036 59 


0.006 194 73 


-0.002 030 70 


-0.005 03312 


-0.005 42743 


C4 


-0.000 610 67 


-0.00147194 


-0.001278 08 


-0.007099 02 


-0.002 00175 


C5 


-0.000 388 13 


-0.003 775 51 


-0.000 934 77 


-0.001 103 32 


0.000100 74 



£$/(Z 2 m/M) 



2 


0.765 698 46 


0.538 869 48 


0.545 667 88 


0.542 471 90 


0.517147 94 


3 


0.840 987 69 


0.562 509 01 


0.569 810 61 


0.582 725 90 


0.524 732 09 


4 


0.879 755 41 


0.576 178 76 


0.582 830 44 


0.608 268 18 


0.529 299 61 


5 


0.903 348 97 


0.584 990 29 


0.590 909 14 


0.625 203 60 


0.532 349 14 


6 


0.919 210 60 


0.591 123 83 


0.596 399 30 


0.637117 06 


0.534 516 27 


7 


0.930 603 34 


0.595 633 60 


0.600 370 03 


0.645 913 88 


0.536130 48 


8 


0.939 18163 


0.599 086 96 


0.603 374 29 


0.652 66102 


0.537 37739 


9 


0.945 873 29 


0.601815 29 


0.605 726 19 


0.657 993 92 


0.538 368 63 


10 


0.951238 85 


0.604 024 79 


0.607 61717 


0.662 312 12 


0.539175 09 


11 


0.955 636 86 


0.605 850 39 


0.609 170 53 


0.665 878 55 


0.539 843 79 


12 


0.959 307 33 


0.607 384 02 


0.610 469 21 


0.668 872 98 


0.540 40711 



1/Z expansion coefficients 



CO 


1 


5/8 


5/8 


5/8+2 9 /3 8 


5/8-2 9 /3 8 


Cl 


-0.491 706 5 


-0.219 6812 


-0.176 992 4 


-0.422 232 8 


-0.082 563 1 


C2 


0.039 6517 


0.100 3376 


0.030 684 4 


0.133 727 


0.045 923 6 


C3 


0.012 972 5 


-0.009 874 4 


0.009 154 5 


0.162 310 7 


0.009 983 7 


C 4 


-0.000 022 6 


-0.006 744 6 


0.003 896 3 


-0.000 445 3 


-0.007945 4 


C5 


0.003 0374 


0.008 9491 


0.001835 8 


-0.099 048 5 


-0.011326 2 



£^ 2 /(Z 2 m 2 /M 2 ) 



2 


-0.923 067 75 


-0.575 065 28 


-0.561902 54 


-0.596 056 62 


-0.552 238 02 


3 


-1.015 027 87 


-0.623 263 22 


-0.591803 40 


-0.687 746 68 


-0.574156 42 


4 


-1.061 762 41 


-0.651489 97 


-0.607 626 64 


-0.742 385 92 


-0.583 335 75 


5 


-1.090 05134 


-0.669 68217 


-0.617 364 45 


-0.776 892 85 


-0.588 425 34 


6 


-1.10901461 


-0.682 312 57 


-0.623 95164 


-0.800 354 77 


-0.591687 77 


7 


-1.122 610 79 


-0.691572 88 


-0.628 70172 


-0.81726102 


-0.593 96708 


8 


-1.132 835 81 


-0.698 645 59 


-0.632 288 26 


-0.829 994 71 


-0.595 653 31 


9 


-1.140 805 12 


-0.704 220 80 


-0.635 09175 


-0.839 919 81 


-0.596 952 97 


10 


-1.147190 96 


-0.708 727 06 


-0.637 343 21 


-0.847868 27 


-0.59798612 


11 


-1.152 422 61 


-0.712 44412 


-0.639 190 98 


-0.854 374 66 


-0.598 827 53 


12 


-1.156 787 03 


-0.715 56217 


-0.640 734 67 


-0.859 79749 


-0.599 526 27 
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TABLE II: The leading relativistic corrections £!£' and £$ for helium-like atoms and their 1/Z-expansion coefficients. The 
analytical results for the coefficient ci for £$ were taken from Ref. [29] for the l 1 ^, 2 3 5', 2 3 Po, and 2 3 P2 states. For the other 
states, this coefficient was evaluated numerically to high accuracy in this work by the same method as in Ref. [29]. The co 
coefficient of for the S states originates from the one-electron recoil effect and is well known from the hydrogen theory. 
For the P states, it contains also the two-electron contribution, which was derived in Ref. [36]. The remaining 1/Z-expansion 



coefficients were obtained by fitting the numerical data for £^ and £^' ■ Atomic units are used 



?(4) 



z 


1 L S 


2 1 S' 


2 3 S 


2 1 P 1 


2 3 P 


2 3 Pi 


2 3 P 2 




/Z 4 














2 


-0.121984 67 


-0.127135 46 


-0.135 279 87 


-0.127 50160 


-0.118 042 52 


-0.123 316 23 


-0.123 729 58 


3 


-0.145 794 73 


-0.13188123 


-0.142 840 30 


-0.130 953 91 


-0.121 128 29 


-0.126 60186 


-0.124 410 72 


4 


-0.163 263 53 


-0.136 297 75 


-0.147 356 68 


-0.133 229 53 


-0.126 667 74 


-0.130 512 05 


-0.125 569 23 


5 


-0.176 048 64 


-0.139 882 98 


-0.150 310 08 


-0.134 786 60 


-0.131522 71 


-0.133 708 51 


-0.126 55782 


6 


-0.185 67419 


-0.142 73756 


-0.152 383 24 


-0.135 91790 


-0.135 444 97 


-0.136 216 54 


-0.127 342 75 


7 


-0.193138 21 


-0.145 028 97 


-0.153 91610 


-0.136 778 46 


-0.138 596 36 


-0.138 199 25 


-0.127 965 82 


8 


-0.199 077 87 


-0.146 895 44 


-0.155 094 63 


-0.137 455 97 


-0.141 156 38 


-0.139 793 29 


-0.128 46748 


9 


-0.203 909 05 


-0.148 439 23 


-0.156 028 56 


-0.138 003 71 


-0.143 26617 


-0.141097 57 


-0.128 87810 


10 


-0.207911 70 


-0.149 734 47 


-0.156 786 71 


-0.138 455 96 


-0.145 029 86 


-0.142 18211 


-0.12921951 


11 


-0.211280 06 


-0.150 835 16 


-0.157414 35 


-0.138 835 82 


-0.146 523 58 


-0.143 096 91 


-0.129 50741 


12 


-0.214 152 65 


-0.15178123 


-0.157 942 44 


-0.139 159 46 


-0.147 803 54 


-0.143 878 26 


-0.129 753 21 


1/Z expansion coefficients 












CO 


-1/4 


-21/128 


-21/128 


-55/384 


-21/128 


-59/384 


-17/128 


Cl 


0.480 139 61 


0.169 478 18 


0.076 935 23 


0.055 403 03 


0.219 768 22 


0.130 428 76 


0.040 638 72 


C2 


-0.636 506 86 


-0.281858 62 


-0.042 775 47 


-0.090 63215 


-0.303 523 35 


-0.162 12941 


-0.047 315 68 


C3 


0.456 314 23 


0.202 919 21 


0.010 473 95 


0.156 412 39 


0.091 746 25 


0.042 468 90 


0.002 244 38 


c 4 


-0.171 179 61 


-0.042 542 10 


-0.004 460 83 


-0.178 042 53 


-0.008 844 33 


-0.004 319 44 


-0.000 236 51 


C5 


0.018 58749 


0.018 86171 


-0.001566 73 


0.059 068 31 


0.015 552 82 


0.007 698 46 


0.003 69105 



£$/(Z 4 m/M) 



2 


-0.134 960 7 


-0.004 3516 


0.005 574 1 


-0.003 655 3 


0.015 596 8 


0.016 6771 


0.012 760 7 


3 


-0.123 759 2 


-0.0016161 


0.011426 9 


-0.008 574 4 


0.026 148 2 


0.026 855 2 


0.019 248 5 


4 


-0.107 6271 


0.002 303 9 


0.015 288 3 


-0.012 139 6 


0.032 666 5 


0.032 185 5 


0.021650 8 


5 


-0.093 784 1 


0.005 792 4 


0.017 933 4 


-0.014 4516 


0.037 580 2 


0.035 768 1 


0.022 926 5 


6 


-0.082 625 7 


0.008 672 2 


0.019 840 8 


-0.015 970 4 


0.0414077 


0.038 388 7 


0.023 736 9 


7 


-0.073 647 


0.011026 9 


0.021276 3 


-0.017 0051 


0.044 450 5 


0.040 395 3 


0.024 304 7 


8 


-0.066 337 


0.012 965 9 


0.022 393 8 


-0.017 736 6 


0.046 915 2 


0.041 981 4 


0.024 727 6 


9 


-0.060 2991 


0.014 581 1 


0.023 287 6 


-0.018 2710 


0.048 946 2 


0.043 266 3 


0.025 0561 


10 


-0.055 2416 


0.015 943 


0.024 018 4 


-0.018 672 8 


0.050 645 5 


0.044 328 


0.025 319 3 


11 


-0.050 950 5 


0.017104 5 


0.024 626 8 


-0.018 982 2 


0.052 086 4 


0.045 219 7 


0.025 5351 


12 


-0.047 267 8 


0.018 105 5 


0.025 141 1 


-0.019 225 6 


0.053 322 8 


0.045 979 1 


0.025 715 5 


l/Z expansion coefficients 












CO 





1/32 


1/32 


-0.020 744 7 


0.069 205 9 


0.055 392 


0.027 764 2 


Cl 


-0.645 040 2 


-0.182 643 4 


-0.078 412 4 


0.002 583 


-0.217113 6 


-0.125 397 2 


-0.025 605 


C2 


0.972 372 8 


0.314 800 3 


0.062 834 3 


0.220104 6 


0.332 302 1 


0.157 607 


0.015 453 2 


C3 


-0.460 0919 


-0.188 443 6 


-0.018 866 9 


-0.3874951 


-0.152 038 8 


-0.096 226 5 


-0.036 663 9 


C4 


-0.040 368 


-0.048 228 2 


0.000 413 8 


-0.046 282 4 


-0.252 878 4 


-0.030 160 4 


-0.017 642 5 
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TABLE III: Bethe logarithm for helium-like atoms with the infinite nuclear mass, hi(ko/Z 2 ), and its first-order perturbation by 
the mass polarization operator, ln(fco)M. Coefficients of the 1/Z expansion of \n(ko/Z 2 ) are also presented. The leading term 
Co is known with a high accuracy from the hydrogen theory. The higher-order coefficients are obtained by fitting the numerical 
data. 

_Z_ f\5 2\S* 2^5 2\P 2^P Ref. 

Hko/z 2 ) 



2 


2.983 865 8618(1) 


2.980118 3651 (1) 


2.977 742 459 29 (2) 


2.983 803 382 4(1) 


2.983 691003 3(2) 






2.983 865 860 9(1) 


2.980118 364 8(1) 


2.977 742 459 2(1) 


2.983 803 377(1) 


2.983 690 995(1) 


[25] 




2.983 865 857(3) 


2.980118 36(7) 


9 Q77 7 A 9 AC\ (\~ \ 
Z.\3 ( I /4z 40 yi) 


z.yoo oUo 40 yo ) 


9 OQQ ftOfl QA f9\ 

z.yoo OyU o4 yZ ) 


ro/ii 

[24J 


3 


2.982 624 563 0(2) 


2.976 363 063 0(2) 


9 Q7^ 8^1 70Q Q9 ( A\ 


9 QK^ 1 01 ^ (\ (9\ 

z.yoo loo uio o yz) 


9 Q°.9 Q^R 7Q8 9 (9\ 

z.yoz yoo / yo z yz) 


4 


2.982 503 0991 (3) 


2.973 976 9112(3) 


9 Q71 7*3^ ^78 QO (7\ 


9 G89 fiQ°. 91^8 {A} 

z.yoz oyo zio o y^) 


9 Q89 A flS 

z.yoz 440 oyo 4 ^o^ 




5 


2.982 591 376 1 (4) 


2.972 388 098 8 (4) 


9 G70 A.9A Qft/1 QO (Sl\ 

z.y < u 4Z4 yo4 yu yo i 


9 Q89 ^/10 1 1 A Q 

z.yoz 04U 114 y ^oj 


9 Q&9 OSQ £0/1 Q f/H 

z.yoz uoy OU4 y ^4^ 




6 


2.982 716 948(1) 


2.971 266 246 4(5) 


9 QfiQ ^7 071 Q (1\ 

z.yoy oo / u / 1 y yo) 


9 QS39 079 71 Q f 9\ 
Z.yoZ U / Z ( Lif yz) 


9 QS1 Q*}8 ^ 

z.yoi ooo yoo o ^oj 






2.982 716 948 (4) 


2.971266 24(4) 


2.969 537 07(1) 


2.982 072 76 (2) 


2.981835 92(3) 


[241 

L J 


7 


2.982 839 085 (3) 


2.970 435 367(1) 


2.968 896 814(1) 


2.981867337(7) 


2.981646 451 (2) 




8 


2.982 948 318 (4) 


2.969 796 528 (3) 


2.968 413 645 (2) 


2.981705 33(1) 


2.981499 939 (4) 




9 


2.983 043 667(8) 


2.969 290 586 (5) 


2.968 036 227(5) 


2.981574 56 (3) 


2.981383 443(7) 




10 


2.983 126 46 (2) 


2.968 880 24(1) 


2.967 733 341 (9) 


2.981466 92 (5) 


2.981288 68(1) 




11 


2.983 198 50 (3) 


2.968 540 85 (2) 


2.967 484 93 (2) 


2.981376 9(1) 


2.98121012 (2) 




12 


2.983 26147(5) 


2.968 255 57(4) 


2.967 277 54(3) 


2.981300 4(2) 


2.981 143 96 (3) 




Coefficients of the 1/Z expansion: 










CO 


2.984 128 56 


2.964 977 59 


9 Qfi/I Q77 

z.yo4 y j i oy 


9 QRO 17fl A7 
z.you O ( O 4 i 


9 Q80 ^7fi A 7 
Z.yoU 0(04/ 




Cl 


-0.012 299 28 


0.040 788 09 


O 097 7^Q A*\ 

u.uz* /oy 4o 


019 OO*} 
U.U1Z UUO oo 


O OOQ «97 Q7 

u.uuy oz i y / 




C2 


0.022 449 74 


-0.016 439 35 


O OOl /19*3 OA 

— u.uui 4zo yo 


OOIO Q89 OK 

— u.uiu yoz uo 


O OO/I 81 O fiO 
— U.UU4 OlU OU 




C3 


0.003 586 19 


-0.012 355 03 


o ooc; q^s K(\ 
— u.uuo yoo do 


O OOO /I £9 1 Q 

— u.uuu 4oz iy 


O 009 A ^7 88 
— U.UUZ 40 i oo 




C4 


-0.002 503 70 


0.005 813 30 


O OOO 1 1 Q 

u.uuu iiy yo 


O OO^ 7«A ^9 
U.UUO I 04 OZ 


O OOO 9^R 70 
— U.UUU ZOO /U 




ln(fc 


M /(m/M) 












2 


0.094 389 4(1) 


0.017 734 4(1) 


0.004 785 54 (1) 


—0.003 553 4 (2) 


0.008 709 5 (1) 






0.094 38(1) 


0.017 734(1) 


0.004 784(3) 


-0.003 538 (6) 


0.008 701 (4) 


[24] 


3 


0.109 539 7(1) 


0.034 210 3(1) 


0.007 852 51 (1) 


-0.006 602 3(2) 


0.0163283(1) 




4 


0.1169197(1) 


0.044 876 8(1) 


0.009 616 61(1) 


-0.007 9512(2) 


0.020199 2(1) 




5 


0.121304 5(2) 


0.052 012 4(2) 


0.010 754 20(1) 


-0.008 629 5(2) 


0.022 479 0(1) 




6 


0.124 212 9(3) 


0.057 053 0(3) 


0.011547 92(1) 


-0.009 015 3(2) 


0.023 973 6(1) 






0.12421 (1) 


0.057 051 (1) 


0.011541(1) 


-0.008 98(1) 


0.023 98(1) 


[24] 


7 


0.126 2831 (4) 


0.060 783 (9) 


0.012 133 20(1) 


-0.009 255 6(2) 


0.025 0273(1) 




8 


0.127 833 6(5) 


0.063 647 0(7) 


0.012 582 68(1) 


-0.009 416 0(2) 


0.025 809 5(2) 




9 


0.129 0371(2) 


0.065 912(1) 


0.012 938 76(1) 


-0.009 528 7(2) 


0.026 412 8(3) 




10 


0.129 998 9(2) 


0.067 746 (1) 


0.013 227 87(1) 


-0.009 6112(2) 


0.026 892 2(4) 




11 


0.130 7851(2) 


0.069 261 (1) 


0.013 467 30(1) 


-0.009 673 6(2) 


0.027 282 2(4) 




12 


0.131439 7(2) 


0.070 534(2) 


0.013 668 86(1) 


-0.009 722 2 (3) 


0.027 605 7(4) 
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TABLE IV: The leading logarithmic QED corrections S^\\og) and ^(log). For the non-recoil correction, we present the 
coefficients of the 1/Z expansion obtained by fitting the numerical data (except for Co which is known analytically). The recoil 
correction is very small for ions with Z > 12, so its 1/Z expansion was not studied. Atomic units are used. 

Z 1 X S 2 1 S 2 3 S 2 1 P 2 3 P 



££>Qog)/[Z* ln(Za)- 2 } 



2 


0. 


.587 967 740 


0.435 225 697 


0.440 118 361 


0.424 690 417 


0.419 620 202 


3 





661366 949 


0.444 453 845 


0.450 745 298 


0.425 056 974 


0.418 277 856 


4 





,702 709 966 


0.450 600 513 


0.456 788 812 


0.425 087 644 


0.418 654 616 


5 


0. 


.729 153 414 


0.454 883 655 


0.460 628 756 


0.425 031723 


0.419 253 068 


6 





.747 506 286 


0.458 013 834 


0.463 274 004 


0.424 963 548 


0.419 808 283 


7 





,760 984 320 


0.460 392 952 


0.465 204 144 


0.424 901229 


0.420 280 937 


8 





,771300 325 


0.462 259 020 


0.466 673 600 


0.424 848 038 


0.420 676 712 


9 





,779 449 465 


0.463 760 328 


0.467 829 292 


0.424 803 420 


0.421008 829 


10 





,786 049 205 


0.464 993 553 


0.468 761 819 


0.424 766 013 


0.421289 730 


11 


0. 


,791502 875 


0.466 024 214 


0.469 530 027 


0.424 734 472 


0.421529 561 


12 





,796 085 045 


0.466 898 202 


0.470 173 774 


0.424 707 664 


0.421 736 263 



1/Z expansion coefficients 



CO 


8/(3tt) 


3/(2tt) 


3/(277) 


4/(3*0 


4/(3*0 


Cl 


-0.659 550 48 


-0.137 744 61 


-0.089 756 44 


0.003 158 46 


-0.036 478 76 


C2 


0.330 586 16 


0.136 452 03 


0.026 693 63 


0.009117 22 


0.051362 86 


C3 


-0.132 768 16 


-0.061297 50 


0.005 48816 


-0.06123169 


0.011 13217 


C4 


0.048 550 42 


-0.000 485 58 


0.001979 36 


0.070 367 56 


-0.001924 41 


CB 


-0.005 834 76 


-0.000 419 68 


0.001082 09 


0.002 994 02 


-0.014 225 78 



£<?(log)/[(m/M) Z 4 ln(Za)- 2 ] 



2 


-1.490 787 8 


-1.087 827 3 


-1.099 5716 


-1.048 033 4 


-1.072 9421 


3 


-1.503 000 7 


-0.999 820 


-1.013 525 6 


-0.931087 7 


-0.972 9291 


4 


-1.4148151 


-0.900 768 2 


-0.913 023 5 


-0.820 369 5 


-0.869 549 4 


5 


-1.281038 5 


-0.795 295 8 


-0.805 624 3 


-0.712 358 9 


-0.765 222 3 


6 


-1.122 852 5 


-0.685 953 6 


-0.694 495 7 


-0.605 509 9 


-0.660 414 8 


7 


-0.950 101 1 


-0.574135 7 


-0.581 1375 


-0.499 185 8 


-0.555 305 8 


8 


-0.767 968 8 


-0.460 648 8 


-0.466 344 1 


-0.393 108 1 


-0.449 990 7 


9 


-0.579 443 


-0.345 987 


-0.350 573 6 


-0.287146 9 


-0.344 527 5 


10 


-0.386 365 2 


-0.230 468 4 


-0.234108 4 


-0.181238 4 


-0.238 954 2 


11 


-0.189 9321 


-0.1143064 


-0.1171318 


-0.075 350 3 


-0.133 296 7 


12 


0.009 045 2 


0.002 350 7 


0.000 232 1 


0.030 5341 


-0.027 573 6 
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TABLE V: The leading nonlogarithmic QED corrections £<^(nlog) and fj^(nlog). For the non-recoil part, we present the 
coefficients of the 1/Z expansion. The coefficient Co is known with a very good accuracy from the hydrogen theory. The 
remaining coefficients were obtained by numerical fitting. The radiative recoil correction is very small for ions with Z > 12, so 
its l/Z expansion was not studied. Atomic units are used. 

Z l x S 2 1 S' 2 3 S 2 1 P 1 2 3 P 2 3 Pi 2 3 P 2 



fL 5) (nlog)/Z 4 



2 


-1.390 282 4 


-1.021756 7 


-1.032 719 5 


-0.999 104 6 


-0.986 487 5 


-0.987 824 8 


-0.987 273 6 


3 


-1.552 423 4 


-1.040 733 9 


-1.055 874 3 


-1.000 028 8 


-0.984 014 5 


-0.985 285 6 


-0.983 645 7 


4 


-1.645 8291 


-1.053 494 3 


-1.068 9514 


-0.999 838 9 


-0.985 557 5 


-0.986 317 8 


-0.983 817 7 


5 


-1.706 668 4 


-1.062 5513 


-1.077 230 3 


-0.999 427 5 


-0.987 480 3 


-0.987 748 5 


-0.984 6170 


6 


-1.749 464 5 


-1.069 274 8 


-1.082 920 6 


-0.999 028 5 


-0.989 184 7 


-0.989 040 9 


-0.985 4376 


7 


-1.7812119 


-1.074 448 4 


-1.087 066 2 


-0.998 687 5 


-0.990 609 4 


-0.990 128 4 


-0.986 162 1 


8 


-1.805 7010 


-1.078 546 


-1.090 218 7 


-0.998 404 7 


-0.9917910 


-0.991033 


-0.986 779 7 


9 


-1.825165 8 


-1.081868 4 


-1.092 695 9 


-0.998 170 7 


-0.992 776 8 


-0.991 788 8 


-0.987303 6 


10 


-1.841008 7 


-1.0846148 


-1.094 693 4 


-0.997 975 9 


-0.993 6076 


-0.992 426 


-0.987 749 6 


11 


-1.854154 5 


-1.086 922 2 


-1.096 338 


-0.997 812 2 


-0.994 315 


-0.992 968 8 


-0.988132 3 


12 


-1.865 237 8 


-1.088 887 5 


-1.097 7156 


-0.997 673 3 


-0.994 923 5 


-0.993 435 8 


-0.988 463 2 



.-(5), 



\/Z expansion coefficients 

c -1.995 417 -1.113 2781 -1.113 2781 -0.996 116 -1.002 747 5 -0.999 4318 -0.992 800 3 

ci 1.658 816 0.325 517 0.191 147 5 -0.017 559 8 0.106 210 2 0.081 192 9 0.059 452 5 

c 2 -1.226 2714 -0.420 456 -0.051603 3 -0.031853 8 -0.146 537 2 -0.108 500 3 -0.085 893 4 

c 3 0.825 843 5 0.327 325 -0.013194 5 0.247 256 7 -0.016 628 5 -0.029 527 6 -0.036 773 9 

c 4 -0.373 062 3 -0.114 627 2 -0.008 9491 -0.334 953 5 0.006 846 2 0.008 2013 0.009 4312 

c 5 0.088 160 0.037 575 2 0.010 984 6 0.085 299 0.018 205 0.016 1412 0.015 502 4 

nlog)/[(m/M) Z 5 ] 

2 3.292 520 2.455 583 2.489 805 2.385 181 2.393 644 2.393 984 2.393 698 

3 2.816 621 1.914 328 1.949 969 1.811294 1.817 918 1.818 182 1.817 751 

4 2.527 383 1.642 421 1.673 220 1.526 466 1.531293 1.531426 1.530 930 

5 2.334128 1.478 228 1.504 411 1.356 337 1.359 930 1.359 969 1.359 460 

6 2.196 193 1.368117 1.390 605 1.243 181 1.245 918 1.245 898 1.245 400 

7 2.092 908 1.289 066 1.308 659 1.162 444 1.164 573 1.164 516 1.164 039 

8 2.012 719 1.229 525 1.246 828 1.101918 1.103 606 1.103 526 1.103 074 

9 1.948 679 1.183 048 1.198 512 1.054 847 1.056 206 1.056113 1.055 685 

10 1.896 366 1.145 753 1.159 714 1.017186 1.018 295 1.018194 1.017 790 

11 1.852 835 1.115 160 1.127 873 0.986 367 0.987 283 0.987177 0.986 794 

12 1.816 049 1.089 607 1.101272 0.960 678 0.961441 0.961334 0.960 971 
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TABLE VI: The ma 6 corrections ^(log) and ^'(nlog) and their l/Z expansion coefficients. Atomic units are used. 



z 


l 1 ^ 


2 1 S 


2 3 S' 


2 1 Pi 


2 3 Pn 
" J u 


2 3 R 


2 3 p 2 


c(e) 


(log)/[Z 3 ln(Za)- 2 } 














2 


0.020 880 865 


0.001698116 





0.000 144 350 











3 


0.031050 719 


0.003 738 928 





0.000 572 351 











4 


0.037 377 475 


0.005 250 650 





0.001013 004 











5 


0.041625 375 


0.006 344104 





0.001386 910 











6 


0.044 658 932 


0.007157 223 





0.001692 317 











7 


0.046 928 852 


0.007 781 230 





0.001 941 610 











8 


0.048 689 367 


U.UUo 1(6 Dlo 


u 


U.UUZ 14 ( (Jyz 


u 


u 


U 


9 


0.050 093 825 


r\ f\i~\0 Of? C 

0.008 671 365 





0.002 318 555 











10 


0.051239 918 


r\ t~\r\c> r\r\r\ noo 

0.008 999 033 





0.002 463 398 











11 


0.052 192 721 


r\ r\r\r\ ri to a ti 

0.009 273 471 





0.002 587 159 











12 


0.052 997 208 


r\ r\r\r\ f- r\r t f T/~\ 

0.009 506 579 





0.002 694 007 











l/Z expansion coefficients 
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Cl 
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n 

u 
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u 


n 
u 


C2 


0.091963 7 


n rns 71 1 4 

U.UOO f 111 


n 
u 


u.uo* yyo t 


n 
u 


u 


n 

u 


C3 


-0.033 400 9 


-fl 014 983 9 


n 


— fl f!39 QQ^ 1 
u.uoi yyo i 


n 

u 


n 

u 




C4 


0.004 804 4 


U.UUO uuu u 




U.UUO OUJ o 


n 
u 


n 
u 


n 

u 


C5 


-0.000 487 4 


-A nnfi 979 (S 

U.UUU ^ 1 ^ u 


u 


n nn9 rp,7 7 

U.UUZ OU 1 1 


n 


n 


n 

u 


c(6) 


(nlog)/Z 6 














2 


2.181233 3(1) 


1.528 981 (2) 


1.536 5931 (1) 


1.489 195 (1) 


1.459 456(1) 


1.460 802 (1) 


1.466 251 (1) 


3 


1.582 4717(1) 


1.016 337(1) 


1.020 440 5(1) 


0.977 073(1) 


0.945 154(1) 


0.937 827(1) 


0.950 303(1) 


4 


1.224 4519 (8) 


0.753 467(1) 


0.754 660 8 (1) 


0.723 937 (1) 


0.691625(1) 


0.676 327(1) 


0.696 280 (1) 


5 


0.989 185 (1) 


0.592 746 (1) 


0.592 238 5(1) 


0.574 526(1) 


0.539 381(1) 


0.517140(1) 


0.544 513(1) 


6 


0.823 3481 (1) 


0.484 094(1) 


0.482 626 2(4) 


0.477 044(1) 


0.437 493(1) 


0.408 810 (1) 


0.443 440(1) 


7 


0.700 317 8 (1) 


0.405 658 (1) 


0.403 650 3(1) 


0.409 248 (1) 


0.364 415 (1) 


0.329 514 (1) 


0.371 247 (1) 


8 


0.605 469 6(1) 


U.o4D o41 


U.o44 UoD Z 


U.oDU UUo [l) 


U.oUy oyy [L) 


fl OAS Q7/1 /I ^ 

U.zDo o ( 4 [L) 
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0.530 1381 (1) 


u.zyy oao (ij 


U.zy ( 430 U (l) 


U.3Z3 13D (I) 


U.ZDD 4oo (I) 


U.ziy 3t>3 (l) 


n 07/1 no 1 1 ^ 
U.A 1 4 y3o (1 ) 


10 


0.468 8716 (1) 


0.262 537 (1) 


0.260 008 6 (1) 


0.294 919 (1) 


0.232 024 (1) 


0.178 826 (1) 


0.241 202 (1) 


11 


0.418 072 5(6) 


0.231 831 (1) 


r\ c\c\r\ o,ot -t /i \ 

0.229 287 1 (1) 


0.272 996 (1) 


0.203 772 (2) 


r\ 1 A A A OO f A\ 

0.144 482 (4) 


0.213 589 (1) 


12 


0.375 2721 (1) 


0.206 144 (1) 


0.203 617 (4) 


0.255 791 (1) 


0.180 178 (2) 


f\ -I 1 A TOO / A \ 

0.114 783 (4) 


0.190 568 (1) 


l/Z expansion coefficients 
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729/114688 





-729/114688 





CO 


-1/8 


-85/1024 


-85/1024 


-0.0792398 


-85/1024 


-0.0672446 


-65/1024 


Cl 


6.342 8979 


3.549 6121 


3.487 5483 


3.099 80 


3.204 5132 


3.11905 


3.059 7402 


C2 


-4.2619 


-1.042 6 


-0.590 


0.068 6 


-0.623 6 


-0.260 9 


-0.1578 


C3 


2.4241 


1.087 


0.155 7 


-0.039 


0.833 4 


0.328 9 


0.2714 


C4 


-2.408 8 


-0.815 7 


0.053 9 


-0.4013 


-0.210 5 


0.0717 


0.1010 
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TABLE VII: The finite nuclear size correction Ef s (the used values of the root-mean-square nuclear charge radii are listed 
in Table VIII) and the higher-order correction £'^ 7+ - ) = + f^ad ■ Contributions to the ionization energy are presented. 

Numerical values of the finite nuclear size correction are scaled by the same factor as for the higher-order correction, in order 
to simplify the comparison between them. 



z 


l 1 ^ 


2 1 S 


2 3 S' 


2 1 Pi 


2 3 Pn 


2 3 Pi 


2 3 Po 


Eta /[ma 


7 Z 6 ] 














2 


3.41(1) 


0.2346 (5) 


0.3039 (7) 


0.00735 (2) 


-0.0914(3) 


-0.0914(3) 


-0.0914(3) 


3 


4.5(1) 


0.393 (8) 


0.47(1) 


0.0166 (4) 


-0.110(3) 


-0.110(3) 


-0.110(3) 


4 


3.16(3) 


U.OUO J 


u.oo i yoj 
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00188 (11 


—0 00797 (41 
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2 


-8.2(4.1) 


-0.43 (22) 


-0.59 (30) 


0.093(46) 


0.37(18) 


0.35(17) 


0.31 (15) 


3 


-9.5(3.2) 


-0.72 (24) 


-0.89 (30) 


0.063(21) 


0.37(12) 


0.35(12) 


0.31 (10) 


4 


-9.6(2.4) 


-0.83 (21) 


-0.97(24) 


0.055(14) 


0.312 (78) 


0.294 (73) 


0.261 (65) 


5 


-9.4(1.9) 


-0.86(17) 


-0.98 (20) 


0.052(10) 


0.266 (53) 


0.249 (50) 


0.219(44) 


6 


-9.0(1.5) 


-0.87(14) 


-0.96(16) 


0.0510 (85) 


0.230 (38) 


0.214(36) 


0.187(31) 


7 


-8.7(1.2) 


-0.86(12) 


-0.94(13) 


0.0500 (72) 


0.202 (29) 


0.188 (27) 


0.162 (23) 


8 


-8.3(1.0) 


-0.84(11) 


-0.91 (11) 


0.0490 (62) 


0.181(23) 


0.167(21) 


0.144(18) 


9 


-8.00(89) 


-0.825(91) 


-0.886 (98) 


0.0477 (54) 


0.163(18) 


0.150(17) 


0.129(14) 


10 


-7.70(77) 


-0.806 (80) 


-0.859 (85) 


0.0463 (48) 


0.148(15) 


0.136(14) 


0.116(12) 


11 


-7.43(67) 


-0.786 (71) 


-0.834 (75) 


0.0450 (43) 


0.135(13) 


0.124(12) 


0.1066(97) 


12 


-7.18(60) 


-0.768 (63) 


-0.811(67) 


0.0435 (39) 


0.124(11) 


0.1132 (99) 


0.0984(82) 



TABLE VIII: Total theoretical ionization energies of n — 1 and n = 2 states in light helium-like ions, in cm 1 . A is the nuclear 
mass number and P c h is the root-mean-square nuclear charge radius. 



z 


A 


Rch [fm] 


l'S 


2 X S 


2 3 S 


2 


4 


1.676 (3) 


-198 310.665 1 (12) 


-32 033.228 734(63) 


-38 454.694 593 (86) 


3 


7 


2.43 (3) 


-610 078.549 (11) 


-118 704.799 71 (79) 


-134044.256 96 (98) 


4 


9 


2.52(1) 


-1241256.625 (45) 


-260 064.342 7 (38) 


-284 740.785 8(45) 


5 


11 


2.41(3) 


-2 091995.58(13) 


-456 261.994(12) 


-490 434.928(14) 


6 


12 


2.470(2) 


-3162 423.60(32) 


-707 370.691 (31) 


-751 130.806 (34) 


7 


14 


2.558(7) 


-4452 723.93 (66) 


-1013 458.475 (66) 


-1066 874.012 (72) 


8 


16 


2.701 (6) 


-5 963 074.2(1.2) 


-1374 588.68(13) 


-1437 720.48(14) 


9 


19 


2.898 (2) 


-7693 708.5(2.1) 


-1790 837.48 (22) 


-1863 745.75 (24) 


10 


20 


3.005 (2) 


-9 644 843.7(3.5) 


-2 262 278.68 (36) 


-2 345 025.61 (39) 


11 


23 


2.994(2) 


-11816 821.8(5.4) 


-2 789 016.99 (57) 


-2 881669.22 (60) 


12 


24 


3.057(2) 


-14209915.2(8.1) 


-3 371 143.82 (86) 


-3 473 772.54 (90) 


Z 




2 1 P 1 


2 3 P 


2 3 Pi 


2 3 P 2 


2 




-27175.771 929(13) 


-29 222.838110 (54) 


-29 223.826 028 (51) 


-29 223.902 466 (45) 


3 




-108 270.881 302 (69) 


-115 812.954 89 (40) 


-115 818.148 68 (38) 


-115 816.058 08 (34) 


4 




-243 787.567 65 (25) 


-257 876.174 4(14) 


-257887.732 4(14) 


-257 872.840 8(12) 


5 




-434 000.42129 (74) 


-455 041.300 3(37) 


-455 057.499 0(35) 


-455 004.840 0(31) 


6 




-679 022.7617(18) 


-707 232.019 2(81) 


-707244.528 3(76) 


-707108.7312(66) 


7 




-978 929.118 5(39) 


-1014 453.504(15) 


-1014 444.829(14) 


-1014153.830(12) 


8 


-1333 765.988 7(74) 


-1376 741.630 (27) 


-1376 682.831 (25) 


-1376131.273 (22) 


9 


-1743 582.370(13) 


-1794154.597(44) 


-1794 003.382 (41) 


-1793 045.585 (35) 


10 


-2 208 407.741 (22) 


-2 266 761.636(70) 


-2 266 460.989 (65) 


-2 264 903.347(54) 


11 


-2 728 303.815 (35) 


-2 794 653.26(11) 


-2 794130.721 (97) 


-2 791723.350 (82) 


12 


-3 303 295.620 (52) 


-3 377 923.51 (15) 


-3377 091.31 (14) 


-3 373 518.82(12) 
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TABLE IX: Comparison of theoretical and experimental transition energies. Units are MHz for He and Li + and cm 1 for other 
ions. Results by Drake are from 2005 for He [31], from 1994 for Li + [37], and from 1988 for other ions [5]. 

Z This work Drake Experiment Reference 



2 3 P -2 3 S'i transition: 



Z 




ova 7^/1 no/i v/q 


07R 7RA nOO/1 7\ 

zvo i o4 uyy^i / ) 


07R 7P-.A HO/1 «7Q fi/OI \ 

z/D / 04 Uy4.0 / o o{AL ) 


M 1 


3 




546 560 686(32) 


546 560 627 


546 560 683.07(42) 


[37] 


4 




26 864.6114(47) 


26 864.64(3) 


26 864.612 0(4) 


[38] 


5 




35 393.628(14) 


35 393.70(8) 


35 393.627(13) 


[39] 


8 




60 978.85(14) 


60 979.6(5) 


60 978.44(52) 


[40] 


10 




78 263.98(39) 


78 265.9(1.2) 


78 265.0(1.2) 


[40] 


Z J i - 
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3 




546 404 980(31) 


546 404 885 


546 404 978.80(51) 


[37] 


4 




26 853.0534(47) 


26 852.04(3) 


26 853.053 4(3) 


[38] 


5 




35 377.429(14) 


35 377.40(8) 


35 377.424(13) 


[39] 


8 




61 037.65(14) 


61 037.7(5) 


61 037.62(93) 


[40] 


2 3 P 2 - 


-2 3 S X 


transition: 








2 




276 732 186.1(2.9) 


276 732 183(17) 


276 732 186.593 (15) 


[1] 


3 




546 467 655(31) 


546 467 553 


546 467657.21(44) 


[37] 


4 




26 867.9450(47) 


26 867.92(3) 


26 867.948 4(3) 


[38] 


5 




35 430.088(14) 


35 430.02(8) 


35 430.084(9) 


[39] 


8 




61589.21(14) 


61589.0(5) 


61 589.70(53) 


[40] 


10 




80122.3(4) 


80121.6(1.2) 


80121.53(64) 


[41] 


2 1 P 1 - 


-2^0 


transition: 








4 




16 276.775(4) 


16 276.77(3) 


16 276.774(9) 


[42] 


2 3 P r 


-2^0 


transition: 








7 




986.36(7) 


986.6(3) 


986.3180(7) 


[43] 



